Mon. Not. R. Astron. Soc. 000,[7]-?? (2005) Printed 5 Febraary 2008 (MN MbX style file v2.2) 



Galactic Warp in the overdensity of the Canis Major Region 



M. Lopez-Corredoira^* 

^ Instituto de Astrofi'sica de Canarias, C/. Via Ldctea, s/n, E-38200 La Laguna (S/C de Tenerife), Spain 



o 
o 

(N 

< 



> 

m 

(N 

O 
O 

6 



Accepted xxxx xxx xx. Received xxxx xxx xx; in original form xxxx xxx xx 



% 



ABSTRACT 

Bellazzini et al. (2006b) claim that Lopez-Corredoira et al.'s (2002) warp model is totally un- 
able to reproduce the Canis Major structure in the red clump stars. However, slight variations 
in the azimuth of the Lopez-Corredoira et al. (2002) warp model, justified by the uncertainties 
in the parameter as well as the local irregularities with respect to the average model, derive a 
result much closer to the observations of the overdensity south/north. The bump of red clump 
stars with rriK — 13-13.5 around I = 241°, 6 = —8.5° and the depth of the Canis Major 
structure are also explainable in terms of the warp with an appropriate extrapolation of con- 
stant height between galactocentric radii of 1 3 and 1 6 kpc, as observed roughly in the southern 
warp, instead of a monotonically increasing height like the northern warp; and the observed 
velocity distribution of stars cannot exclude the warp possibility. A warp model is therefore 
still a possible explanation of the Canis Major overdensity, and the hypothesis of the existence 
of a dwarf galaxy is unnecessary, although still a possibility too. 
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1 INTRODUCTION 



Bellazzini et al. (2006b, hereafter B06b)^ have claimed that the ex- 
cess of red clump stars in southern with respect to northern galactic 
latitudes between / = 200° and 280° is most probably associated 
with a new dwarf galaxy, namely the Canis Major (CMa) galaxy, so 
that it is far from being explained in terms of the known characteris- 
tics of the Galactic warp model derived from the parameters given 
by Lopez-Corredoira et al. (2002, hereafter L02); they conclude 
that the warp "is totally unable to reproduce the CMa structure" . 
Previous papers (e.g. Martin et al. 2004, Martinez-Delgado et al. 
2005a) have presented other proofs in favour of the existence of the 
dwarf galaxy and against the warp possibility. The discussion on the 
validity of the proofs is a topic of heated debate nowadays, with au- 
thors like Momany et al. (2004, 2006 [hereafter M06]) arguing that 
these proofs can also be reproduced by the Galactic warp+flare. 

Here, I also want to contribute to the debate showing that 
slight variations in L02 warp model derive a result much closer to 
the observations of the overdensity ( i|2.1> . The bump of red clump 
stars with rriK = 13-13.5 around I = 241°, b = -8.5° (§2.2) and 
the depth of the Canis Major structure (§2.3) are also explainable in 
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terms of the warp. I shall make a few comments about the velocity 
distribution of its stars (§2.4). 



2 GALACTIC WARP VS. DWARF GALAXY 
HYPOTHESIS 

2.1 Canis Major overdensity 

The height of the L02 warp over the plane defined by the central 
disc as a function of the galactocentric distance and azimuth is 



„(i?,<^)H = CwRikpc) 



(1) 



with ew = 5.25 ± 0.5, <^w = -5 ± 5° and Cw = 1.2 x 10-^ 
(for the given values of ew = 5.25 and (pw = —5°; the normal- 
ization will change if ew and/or (p^ vary). B06b calculate in their 
fig. 6(down) the south/north overdensity of the red clump stars, 
p{S) — 1.2p(N) due to this L02 warp and the flared disc parameters 
given in L02^. The factor 1.2 was put by B06b to compensate for 
the average south/north asymmetry. The asymmetry between south 
and north is not constant, and it is a very bad approximation to take 
it as B06b does, but I use their expression in order to show that I 
can roughly reproduce the results of their fig. 6 with the warp. We 
confirm their results in Figure^). Effectively, the maximum of the 



There is an error in L02: "+15" pc was written in eq. (20) instead of 
"—15"; the connect value is —15 pc because the Sun is "above" the plane. 
The calculations in L02 were coiTectly canied out with — 15, and the error is 
just typographical. In any case, this does not affect too much to the analysis 
of the warp at high values of R. 
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overdensity is around 270° as they claim (and contrary to tlie claim 
by M06, who used a different maximum definition and different 
constraints on region selection), somewhat far from the centre of 
the observed CMa structure at / fa 244° (B06b, Fig. 6, top). 

L02 only give a warp formula for R ^ 13 kpc; beyond 13 kpc, 
different extrapolations are possible. We also make the same plot 
with a different extrapolation of the L02 model over the range 13 < 
R{kpc) < 16, with a constant height of the warp, 2^,(13 kpc < 
7? < 16 kpc) = Zw{13 kpc), which is more similar to the real 
gas southern warp (see §2.2). This is shown in Fig.0)); the result 
does not change too much. As noted by B06b, the extrapolation 
beyond 13 kpc is not too important because the CMa feature occurs 
approximately at this galactocentric distance. 

However, some attention should be paid to other parameters 
of the warp. L02 have given an approximate model of the "aver- 
age" warp for the whole sky assuming north-south symmetry and 
a power law for the amplitude of this warp. This assumption is 
just a first-order approximation since, as is well known, the Galac- 
tic warp is not symmetric (Burton 1988; Voskes & Burton 2006, 
Levine et al. 2006); our warp is somewhere between an L-warp and 
a S-warp rather than being a pure S-warp. Moreover, the position 
of the Sun also affects our perspective of the northern and south- 
ern warp. Also, there are some other parameters of the disc (galac- 
tocentric distance of the Sun, height of the Sun above the plane, 
flare parameters, scale length, scale height at the Sun galactocen- 
tric distance, etc.) whose variations affect the result. Furthermore, 
the thick-disc component (Cabrera-Lavers et al. 2005), which was 
not included, could also produce some small variations. 

I am not going to explore all the set of parameters here but just 
the azimuth, (fj^, of the warp; that is, the galactocentric angle in 
which the warp has null amplitude. It was checked that the galactic 
longitude of the maximum p{S) — 1.2p{N) is more sensitive to 
the changes in (j)^ than changes in other parameters. For instance, 
changes in ew (with the corresponding change in Cw to preserve 
the amplitude at high R) mainly affect the distance of the maximum 
overdensity rather than its galactic longitude; also, other changes in 
the disc or flare have less effect than changes in 0™ with galactic 
longitude. Figures 0;-e show the density maps for angles </)„ — 
0, +5°, +10° instead of the original value 0^ — — 5 ± 5° given 
by L02. As observed the variation in the position of the maximum 
is notable. With (fj-w ~ +5° {2a from the value given by L02) 
the angle of the maximum is 248°. Just to see that the position 
of the maximum is not so sensitive to slight changes in the shape 
of the warp, we also performed the calculations with Drimmel & 
Spergel's (2001) formula for the warp amplitude (25[R{kpc) — 7]^ 
instead of CwR"" in expression Q; the result was also 248° (Fig. 
Q^. This galactic longitude is not very far from the observed value 
of 240-244° by B06b, and its distance of 5.4 kpc is somewhat less 
than the value r = 7.2 kpc given by B06b. The distance r of the 
maximum could be changed if the flare, warp and scalelength were 
rescaled, or if they showed a different dependence on R from those 
given by L02. For instance, a large scalelength for the disc and the 
flare, and a lower exponent ew places the maximum farther away. 

These plots may be compared in fig. 6 (top) of B06b. The sim- 
ilarity is quite high (although not for r larger than 9-10 kpc since 
in this case the B06b red clump maps are highly contaminated by 
the dwarf population and are very inaccurate). Taking into account 
that we have only modified one parameter, the result is not so bad. 
An angle of (j>nj = +5° is only 2a from the mean value in the 
fit of the average symmetric warp in the star counts. Other models 
of the average warp also give similar values: ^ly = (Freuden- 



reich 1998; Drimmel & Spergel 2001; Robin et al. 2003) or -1-15° 
(Yusifov 2004, M06). 

We must also consider that the error of 5 degrees for (j)w in 
the analysis of L02 is merely statistical, but there might be further 
systematic errors due, for instance, to errors in the assumed lumi- 
nosity function of stars, or deviations from equation Q. As M06 
say, some uncertainties of the L02 model could also come from 
gaps in the region around I = 270° in the data used for the fit. 
However, M06 are not correct when they claim that the L02 warp 
is affected by contamination of dwarfs for red clump stars, because 
L02 do not use red clump stars for the determination of the param- 
eter of the warp but only the total star counts. 

As has been said, given the irregularities of the southern warp, 
it is quite possible that an extra shift in the azimuth in the south- 
ern with respect the average L02-warp might occur. Indeed, the 
analysis by Voskes & Burton (2006) gives precisely the best fit 
for (j)w = +5° (warp maximum at (j) — 95°) for their gas warp. 
The asymmetries in the gas warp (similar to the stellar warp ac- 
cording to L02) are shown by Voskes & Burton (2006) or Levine 
et al. (2006), which should serve to indicate that there is no unique 
large-scale symmetric warp model of the type eq. 0; therefore, we 
cannot extrapolate the exact result at around I ~ 240° to the whole 
Galaxy. Figure 1 1 of Levine et al. (2006) shows how the line of 
nodes defined by tfj^ varies depending on the galactocentric dis- 
tance and from north to south. Indeed, Levine et al. (2006) explain 
the asymmetries with a warp dependence 



z^{R, 0) = Wo{R) + WiiR) sm[(t> - MR)] 
+W2{R)sm[2(t>-MR)]- 



(2) 



The fact that (pi ^ 02 [the difference is up to 12 degrees according 
to Levine et al. (2006)] and depending on R, causes the equivalent 
(pw in expression not to be constant. The first term, the mode of 
m = 1, is dominant for R < 13 kpc, so eq. Q can be a relatively 
good approximation, but if we aim to provide an accurate expla- 
nation for all the deviations from it, perhaps eq. Q is insufficient. 
Here we are not making a fit of the parameters of the global warp 
; instead, we are showing that within the uncertainties the warp is 
compatible with B06b's plots. 

M06 claim that a global and regular warp signature is 
traced to Galactocentric distances of at least ~ 20 kpc, and the 
north/south asymmetry, apart from the fact that (j)w 7^ 0, is due to 
the chance location of the northern warp behind the Norma-Cygnus 
spiral arm. Perhaps this spiral arm is responsible for the asymme- 
try, but data such as those illustrated in fig. 2 of Levine et al. (2006) 
show something else: a lower amplitude for the southern warp for 
-R > 13 kpc (see §2.2), and this is indeed necessary to explain the 
distribution of sources along the line of sight, something which was 
not successfully done by M06. 



2.2 Bump 

Another comment concerning B06b paper is their claim that the 
bump in their fig. 1, in the region 238° < I < 244°, -11° <b < 
— 6° , is due to the dwarf galaxy. In this case, they compared it with 
the predictions of another warp -1- flare model (Robin et al. 2003) 
and failed to reproduce it. However, the bump could be explainable 
in terms of the warp with an appropriate extrapolation for values of 
R > 13 kpc. In §2.1, we saw that the extrapolation is not relevant 
for the position of the maximum overdensity. However, it is rele- 
vant if we want to analyse the sources beyond 13 kpc. Taking the 
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southern warp equal to the northern warp is inappropriate because 
there is an abundance of data for the gas emission of the Galaxy that 
show the asymmetry. At _R > 13 kpc, the m = 2 mode becomes 
important and the extrapolation of the southern and northern warp 
are not equal. Burton (1988) shows that the southern warp is ap- 
proximately of constant height between 1.6i?Q and 2Rq (beyond 
2Rq is unimportant for our analysis) instead of the monotonically 
increasing northern warp. Something similar is observed in Voskes 
& Burton (2006, fig. 16), and we know that the gas warp and stellar 
warp are similar (L02, M06), so the adopted approximate extrapo- 
lation for stars is justified. 

Compare Fig.|2|with fig. 1 of B06b. In Fig.|2| we see a bump 
with a maximum peak at ttik = 13.4 (equivalent to R = 15.4 
kpc), close to the peak at niK = 13-13.5 obtained in B06b. Note, 
however, that only the red clump giants are plotted in Fig.|2| while 
fig. 1 (right panel) of B06b includes all contaminants, especially 
dwarfs over niK = 13.5. Indeed, the decrease in counts beyond 
rriK = 13-13.5 is due to the end of the warp at i? >16 kpc, and, 
of course, the new increase in the counts around rriK = 13.6-13.8 
in fig. 1 of B06b would be due to the dwarf contamination. 



2.3 Colour-magnitude diagram and Canis Major depth 

Figure 7 of B06b gives a comparison of real and synthetic colour- 
magnitude diagrams. This plot is not easy to analyse by eye because 
it includes CMa and foreground stars all together. And it depends 
on the model of the warp used, so we should not say that since one 
model of the warp does not fit, no model of a warp will fit. B06b 
see in their data a similarity with a Gaussian distribution of stars 
with an r.m.s. of 0.8 kpc. Could the warp produce such a colour- 
magnitude diagram? The value of 0.8±0.3 kpc that B06b take is 
from Martinez-Delgado et al. (2005a,b), who claim with their result 
that the warp hypothesis is difficult to reconcile with such a line-of- 
sight depth in the main sequence of a colour-magnitude diagram. 

On the one hand, the Martinez-Delgado et al. (2005a) paper 
contains errors: i) the Gaussian fit of their fig. 3 gives (JMS,totai ~ 
0.57 mag; once the contribution of intrinsic broadening (0.19 mag) 
and latitude dispersion (0.29 mag) are subtracted, it would give 
o'MS,CMa = 0.45 mag (equivalent to 1.6 kpc at d = 7.9 kpc 
instead of 0.8 kpc as they calculate); ii) a Gaussian distribution of 
bins of constant magnitude does not give a Gaussian distribution 
in the density distribution along the line of sight (the second dis- 
tribution is proportional to the first multiplied by a factor 1/d^), 
so we cannot translate a in magnitudes into a in distances, that is, 
a — 1.6 kpc would be the r.m.s. of the function r'^p{r), not the 
r.m.s. of p{r); iii) the assumed constant distribution for the under- 
lying Milky Way stars (dotted line in fig. 3 of Martinez-Delgado et 
al. 2005a) is not even a good first order approximation. 

On the other hand, if we forget the analysis by Martinez- 
Delgado et al. (2005a) but take as correct their data in their fig. 3 
(a FWHM~ 1.6 magnitudes in the distribution in bins of constant 
magnitude), we can compare them with our predictions of the warp 
and see that it is not so different with respect to the predictions of 
the warp (see Fig.|3j: 2.2 mag (2.4 mag if we took into account the 
intrinsic broadening and the latitude dispersion) instead of 1.6 mag. 
The prediction of the L02 warp gives a somewhat broader distribu- 
tion, possibly because the radial dependence is not very accurate. 

M06 claim instead that this structure is a spiral arm. Again, I 
do not agree with the argument by M06 although I agree with their 
general conclusion that the warp can explain the observed facts. No 
spiral arm is needed; it is just a question of a wrong calculation of 



the thickness by Martinez-Delgado et al. (2005a) and a warp with 
appropriate extrapolation over R > 13 kpc. 

The population attributed by Martinez-Delgado et al. (2005a) 
to be a 1-2 Gyr old population in the blue plume of intermediate- 
age open clusters belonging to the CMa dwarf galaxy is indeed 
a young stellar population (^ 100 Myr) of the Galactic spiral 
arms in the background of open clusters, not placed in the puta- 
tive CMa galaxy (Carraro et al. 2005). The reply of B06b to Car- 
raro et al. (2005) seems insufficient. Another paper by Bellazzini 
et al. (2006a) claims that the metallicity of the core of CMa is 
—0.4 $; [M/H] ^ —0.7, a relatively old population; however, this 
is also within the expectation for the outer disc; for an R — 13.1 
kpc, [Fe/H] ~ —0.57 is expected for the Galactic disc according 
to the metallicity distribution by Cameron (1985). 

In conclusion, I do not see in the analysis of colour-magnitude 
diagrams of the Canis Major region any conclusive proof that we 
are observing a population different from that of our own warped 
Galaxy. 



2.4 Velocity of the CMa stars 

The bimodal distribution in the radial velocity of M-stars (Martin et 
al. 2004), presented as a proof that Canis Major is not the warp, re- 
flects two kinds of origins for the sources: one was artificially pro- 
duced by template issues resulting from a fluctuating line spread 
function asymmetry during the different observing nights, as rec- 
ognized by the authors in a later paper (Martin et al. 2005); and 
the other peak can be reproduced by the Galactic rotation (M06). 
In any case, even if M06 were wrong, the kinematics of the warp is 
somewhat complex and unknown, so we cannot discard it. Another 
recent claim of measured motion of CMa perpendicular to the disc 
(Dinescu et al. 2005; Martinez-Delgado et al. 2005b) should not 
be considered as inconsistent with the expected motion of the warp 
because indeed we do not know very much about how the warp 
was formed or its subsequent evolution, whether it is steady or still 
oscillating with respect to the plane — there is no unique scenario, 
and there are at least four possible hypotheses of warp formation 
(Castro-Rodriguez et al. 2002, Sect. 1), each one offering different 
predictions on its motion — or whether the northern and southern 
parts should have similar kinematics (given the north-south asym- 
metry, it is possible that we cannot compare them). Moreover, as 
M06 state, the stars selected to measure the proper motions might 
be contamination not associated with CMa, and "the expected warp 
signature can be and is compatible with negative vertical velocity ". 



3 DISCUSSION 

The L02 model with the modified (j)^ ~ +5° fails slightly to repro- 
duce the CMa feature but not by so much (and we must bear in mind 
that the method of producing the maps also involves certain errors, 
since the red clump stars used are contaminated by dwarfs, late- 
type giants and other spectral types in different ratios depending on 
the line of sight). The bump of red clump stars with rriK = 13- 
13.5 around I = 241°, h = -8.5° or the depth of the Canis Major 
structure are also explainable in terms of the warp (with the ap- 
propriate extrapolation between 13 and 16 kpc of constant height, 
as observed). The blue plume in the colour-magnitude diagram is 
explicable in term of the spiral arm population. The velocity dis- 
tribution of the stars cannot be a proof to exclude it is a warp. The 
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question now arises as to whetiier it is absolutely necessary to in- 
voke the existence of a new dwarf galaxy to explain the red clump 
stars. 

The two options (warp or dwarf galaxy) are usually chosen de- 
pending on the methodology of analysis. Those authors who prefer 
the dwarf galaxy hypothesis assume a fixed model of the warp and 
tend to think that any departure of this model is due to the existence 
of the new galaxy. However, we must always bear in mind that the 
predicted warp features depend on the parameters of the disc, the 
warp itself, the stellar population, the kinematics, etc.; and all this 
knowledge is not so accurate as to allow a perfect agreement with 
all the data, specially for the warp. Those authors who prefer the 
opposite hypothesis claim that whatever you observe is the warp, 
using an ad hoc model of it (as the case in the present paper with a 
modified (j>w)- Perhaps none of the methodologies is appropriate. 
I am neither in favour nor against the dwarf galaxy hypothesis. It 
is quite possible that CMa is a dwarf galaxy, but due to the prox- 
imity of the warp feature, for which we do not have very accurate 
information, it is difficult to disentangle both effects. 
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Figure 1. a) Subtracted density maps (first contour: 1000 star/kpc^ and following contours every 500 star/kpc^; length units kpc; Sun at the center, the cross 
points out the Galactic center, the square indicates the maximum of the density map in B06b data; the filled circle indicates tha maximum of the present 
density map; the dashed circles mark the regions at distances 13 and 16 kpc respectively from the centre of the Galaxy) calculated in the same way as B06b: 
p{S) — 1.2p{N) integrated over 5° < |6| < 15°, excluding \z\ < 0.5 kpc; with the warp model by L02 (Rq = 8 kpc) for R ^ 16 kpc and no warp beyond 
16 kpc. Maximum at I = 267°, distance from the Sun: 10.0 kpc. b) Same as a), but with constant height of the warp 2„ for 13 < R < 16 kpc, and null 
beyond. Maximum at I = 267°, distance from the Sun: 9.5 kpc. c) Same as b), but with ij>w = as azimuth of the warp. Maximum at I = 255°, distance 
from the Sun: 5.7 kpc. d) Same as b), but with (pyy = +5° as azimuth of the warp. Maximum at / = 248°, distance from the Sun: 5.4 kpc. e) Same as b), 
but with (pw = 4-10° as azimuth of the warp. Maximum at I = 248°, distance from the Sun: 5.4 kpc. f) Same as b), but with (pw = +5° as azimuth of the 
warp, and the Drimmel & Spergel (2001) formula for warp amplitude. Maximum at Z = 248°, distance from the Sun: 5.4 kpc. 
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Figure 2. Prediction of the warp model using L02 parameters fertile counts of the Fig. 1 in the paper by B06b (238° ^ / ^ 244°, —11° ^b^— 6°) with 
extrapolation of constant height in 13 < i? < 16 kpc and null beyond. 




Figure 3. Prediction of the warp model using L02's parameters, except for (j>w, which was changed to +5°, for the counts in fig. 3 in the paper by Martinez- 
Delgado et al. (2()05a) (I = 240°, b ^ —8°) with an extrapolation of constant height for 13 < i? < 16 kpc and zero beyond. 



